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A path-integral for the t-J model in two dimensions is constructed based on Dirac quantization, 
with an action found originally by Wiegmann (Phys. Rev. Lett. 60, 821 (1988); Nucl. Phys. B323, 
311 (1989)). Concentrating on the low doping limit, we assume short range antiferromagnetic 
order of the spin degrees of freedom. Going over to a local spin quantization axis of the dopant 
■ fermions, that follows the spin degree of freedom, staggered CP^ fields result and the constraint 

' against double occupancy can be resolved. The staggered CP^ fields are split into slow and fast 

, modes, such that after a gradient expansion, and after integrating out the fast modes and the 

^Sj ' dopant fermions, a CP^ field-theory with a massive gauge field is obtained that describes generically 

incommensurate coplanar magnetic structures, as discussed previously in the context of frustrated 
quantum antiferromagnets. Hence, the possibility of deconfined spinous is opened by doping a 
coUinear antiferromagnet. 
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I. INTRODUCTION 



High temperature superconductivity (HTS) remains an unresolved problem in spite of an enormous research effort 
over more than twenty years [3]. On the theoretical side, a number of phenomcnological theories were proposed 
d, H, 0] that may describe certain aspects of the experimental findings, however, they rely on fundamentally different 
' assumptions, so that it is, at this point, difficult to assess their validity. Also approaches based on microscopic 
^ models like the t-J one were advanced 0], based on the so-called slave boson formulation, where electrons are split 
I into separated spin (spinon) and charge (holons) degrees of freedom. Such a separation introduces a local gauge 
' O ' invariance, that renders mean-field like approximations particularly troublesome. A common feature of the above 
mentioned approaches is the difficulty to connect in a controlled way the states proposed for the doped case to the 
undoped one. 

In fact, the theoretical description of the evolution from a Neel state in the parent compounds towards the state 
of a doped quantum antiferromagnet, lies at the heart of a theoretical understanding of HTS. A phenomenological 
description of doped antiferromagnets was first given by Shraiman and Siggia in a series of seminal papers 0, S ■ 
The picture emerging from a semiclassical treatment of mobile holes in an antiferromagnet corresponds to a coplanar 
twist of the spin background that gives rise to a dipolar field centered on the dopant holes. Similar coplanar structures 
^\ . are also expected in frustrated quantum antiferromagnets [l^, [0, 112] ■ 

The field-theoretic treatment of frustrated quantum antiferromagnets showed on the one- hand, that an 0(4) symme- 
^ ^ try is dynamically generated, and in two dimensions at temperature T = long-range order can set in either through 
' a first order transition or a second order one with exponents corresponding to an 0(4) non- linear a model [Eini. On 
] the other hand, when the effective theory is formulated in terms of CP^ fields, the presence of deconfined spinous can 

■ be induced from the fact that the corresponding gauge field becomes massive [13; [3 13 • Yet, an explicit connection 
^ ■ between doped antiferromagnets and the effective theories for frustrated quantum antiferromagnets is missing. 

, The experimental observation of incommensurate peaks in neutron scattering experiments on La2-a:Sr3:Cu04 

■ (LSCO) [ill and on YBa2Cu307 (YBCO ) [l^ [TtI. 19*1 revived the interest on coplanar magnetic configurations in 
^ , cuprates. A number of theoretical works [20l. l2ll |22. 23, 24] focused on material specific descriptions of the magnetic 

structures at low doping, mostly based on the phenomenology developed by Shraiman and Siggia. Although these 
works were very successful in interpreting the experimental results, demonstrating that the used phenomenology 
describes very well the generic features at low doping, they did not make an explicit connection to a microscopic 
model. 

An alternative way to describe the low energy physics of doped antiferromagnets is the development of an effective 
field-theory based on the symmetries of the system (global SU(2) and U(l) in the case of a doped quantum antifer- 
romagnet) and the possible spontaneously broken symmetry phases. Such a path, reminiscent of chiral perturbation 
theory was followed starting with the Hubbard model, as a representative one for doped antiferromagnets [2^. At 
the moment the predictive power of the derived effective action is not clear. 

The most direct approach for treating a microscopic model like the Hubbard [2^ or t-J [13] models encountered 
until now a number of difficulties. While on the one hand, a mean-field treatment of the t-J model based on the 
slave boson approach 5] led to a qualitative understanding of various experimental results, like the existence of a 
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pseudogap, it is difficult to assess tfieir reliability due to the uncontrolled nature of the mean-field approximation. The 
inclusion of fluctuations of the gauge fields remains until now difficult, since in this case their coupling to the matter 
fields is strong. Although their treatment led to a qualitative description of the doped phase, still progress would be 
desirable on a more quantitative basis [p. On the other hand, numerically exact results could be only obtained by 
diagonalization in rather small clusters |28| so that their interpretation remains inconclusive when contrasted to the 
low energy behavior in HTS. For larger sizes, variational Monte Carlo techniques indicate that the t-J model supports 
superconductivity 29]. Still, in spite of the accuracy of the method, further insight in the system is still missing. 
Unfortunately, large scale quantum Monte Carlo simulations are hindered by the so-called minus sign problem that 
affects the simulation of doped antiferromagnets. 

In view of the described situation, a controlled analytical treatment of a microscopic model is desirable, in particular 
with the possibility of examining the low doping regime, such that the change from an antiferromagnet with long- 
range order at zero doping to a doped situation can be followed in detail. A first step in that direction was made by 
Wiegmann [s^], who obtained an action for the t-J model based on the coherent-states method [3l|. The same action 
was found by using a supersymetric version of the Faddeev-Jackiw symplectic formalism [s^] applied to the t-J model 
[3^, with explicit expressions for the measure of the path integral. 

In the present work we use alternatively the well known procedure of Dirac quantization for constrained systems 
[S^ l to set up the path integral, recovering the results in Ref. [3^, as shown in Sec. [Ill There we shortly review the 
t-J model and its representation in terms of Hubbard X-operators [H, [s^ [s^], that fulfill the graded Lie algebra 
Spl(2,l) |30, 31, 38]. Introducing a series of primary constraints and on the basis of the action proposed by Wiegmann 
[30,|311, all the constraints of the theory are determined. It turns out that only second class constraints appear, such 
that they can be solved by inverting the matrix of constraints, and the Dirac brackets reproduce the algebra of the 
X-operators. On this basis, the path integral can be set up. Passing from the representation by X-fields to real vector 
fields for the spin degrees of freedom and Grassmann variables for the dopant holes [l^, [s^ , we arrive at the action 
that will be the starting point for a long-wavelength expansion. In order to proceed further, we restrict ourselves to 
the low doping limit and neglect terms quadratic in the density of dopant holes. 

The long- wavelength expansion is performed in Sec. lIIIl where a staggered CP^ representation is introduced for the 
spin-fields. Using the CP^ representation it is possible to resolve exactly the constraint against double occupancy [30[, 
that is in general the stumbling block for a controlled treatment of the model. Slow and fast modes of the CP^ fields 
are identified in the same spirit as done for vector fields in quantum antiferromagnets [s^ . The effective action for 
the magnetic degrees of freedom is reached after integrating out the fast CP^ modes and the fermions. The resulting 
effective field-theory corresponds to a CP^ model with a massive gauge field, as was generally discussed in the context 
of frustrated quantum antiferromagnets [isl, [l3]. In the present work, however, we obtained the explicit doping 
dependence of the coupling constants. In Sec. IVIII we discuss the obtained results that open the possibility of having 
deconfined spinous by doping. Some intermediate results are presented in the appendices that may be helpful for 
readers interested in reproducing our results. 



II. DIRAC QUANTIZATION OF THE t-J MODEL 



We introduce first the t-J model and its representation in terms of so-called X-operators that operate only in the 
subspace without doubly occupancy. After discussing shortly the algebra they fulfill, we delineate the procedure of 
Dirac quantization. 



A. The t-J model and X-operators 



The t-J model is defined by the following Hamiltonian in second quantization: 

<i,3> <i,j> ^ ' i 

with = (1 - nj_o.)c|^, hi = Yl,a^i<j^ia^ Si = J2 4a'^<r<r'^ia'^ whcrc the operators cj^ and c^^ denote canonical 

creation and annihilation operators, respectively, for electrons at site i and spin indices a — zt. The operators 
and Cjg. project out doubly occupied states. Then, t gives the hopping amplitude, J is the antiferromagnetic exchange 
coupling and /i the chemical potential. The symbol < i,j > restricts the sums to nearest neighbors. The Hamiltonian 
p]) was obtained from a multiband model [27( and represents the minimal model for cuprates. 
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Introducing so-called X-operators [35|, |36|, |37| defined as 



X, 



a(3 



I a» >< A I , 



with ai — 0,a for site i, the Hamiltonian becomes a bilinear form in such operators. 



H, 



t-j 



^3 



The X-operators fulfill the following graded algebra 



(2) 



(3) 



(4) 



with — (+) corresponding to a commutator (anticommutator). Anticommutation relations appear only when both 
operators are fermionic. Furthermore, we have the completeness condition 



(5) 



A further insight into the graded algebra above can be gained by considering the commutation and anticommutation 
relations of the even (bosonic) and odd (fermionic) parts. Denoting the even generators by Qm, m ~ 1,2,3 and i3, 
and the odd ones by f/i, i = 1, . . . , 4, the commutation relations of the Spl(2,l) algebra are |38| | 



[Qm, Uo.] = l^T'^o^Uf, , [B, [/„] = ie^, 



1/3 

where the 4x4 matrices (T™, C, and e are defined as follows 

C 



cr" 
cr™ 



C 
C 



1 

-1 



(6) 



(7) 



with (7™, m ~ 1,2,3 the Pauli matrices and C = ir^. In the irreducible representation corresponding to the X- 
operators, the generators above look as follows. For the even sector we have 



Qi 

0,3 



X- 
X- 



, Q2 



4(x- 



X- 



B=-iX- 



X- 



(8) 



The first three generators are those of SU(2) while the last one corresponds essentially to particle number. Here and 
in the following we eliminate using the completeness relation ([5]). For the odd sector we have 



Ui = 1+° , U2 = X^" , U3 = X"' , Ui = -X"+ 
As shown in Ref. [ssl ]. the Casimir operator quadratic in the generators is given by 

K2 = Q^-B^ + luCU, 



(9) 



(10) 



that in terms of the X-operators looks as follows. 



■x-+x^ 



X' 



1 2 



'^x^^ - X- 

X'^- X-° + X^ 

Then, it is easily seen, that this Casimir operator has eigenvalue zero in the present irreducible representation. 



1 [ -x+"x"+ - x-"x'- 



(11) 
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B. Dirac quantization for the t-J model 



In the following we consider a classical system with the Lagrangian found by Wiegmann |30l . |31[ , expressed in terms 
of X-fields with, as usual, complex fields corresponding to bosonic operators in Sec. Ill Al and Grassmann fields 
for fermionic ones. 



where the following definitions were introduced 

p.,^X^+Xt^ + XrX^\ u,^X++-Xr-. v,^X+++Xr-- (13) 

Furthermore, we take into account the set of primary constraints found in the supersymetric extension of the symplectic 
formalism introduced by Faddeev and Jackiw 33], where we omit the site indices. 

<^(i) = X++ + X— + p-1 , 

^(2) ^ x+-X-^ + \u'-[l-lvy + p, 



X++ 

^(10) ^ x+"X-+-X-OX++ . (14) 

The constraints are imp osed by setting 0^'^) = and the order of the labels a is such that the supermatrix of constraints 
has a normal form [40l. [4l| . The special choice of 0^^^ was made in order to have a simple expression for the measure 
of the path integral obtained at the end. The constraint (j)^^'' (together with (Z)*-^^) corresponds to the completeness 
relation ([5]), while (p^^^ comes from the Casimir operator (jlip having eigenvalue zero. Furthermore, i/)'-^-' relates the 
empty sites with fermionic holes, and the rest of the constraints above define new product rules instead of the ones 
obeyed by the operators in Sec. Ill Al [33*1. 

There are other nine primary constraints resulting from considering the canonical momenta 

n"^ = , (15) 

where drL/dX"^ refers to the right derivative, the X"^ and its derivatives being elements of the Berezin algebra [i^l 
(see Appendix [X| . The additional nine constraints are listed in Appendix [XI 
Once the constraints are determined, a Hamiltonian 

i7 = ^ n^'^X"'' - L (16) 

can be in principle obtained, where we have to take into account the constraints. This can be done introducing 
Lagrange multipliers A that enter in the equations of motion for any observable f{X, H). In particular, if the primary 
constraints are required to apply at any time, we should require that 

0^"^ = {H,.^("'}+Ab{#\(/.(°)} = 0, (17) 

where the curly brackets denote now Poisson brackets and summation over repeated indices is assumed. It is under- 
stood that the constraints are applied after the derivatives in the Poisson brackets are calculated. For the Poisson 
brackets among the constraints we have the general form ^] 



E 



dxf dnf dxf duf 



(18) 



where P^(a) is the parity of 4'^'^\ and deip'-'''' / dll'^^ refers to the left derivative [i^l, that is defined in a similar way 
as the right derivative in (fT51) (see Appendix [X| . If the matrix of constraints } is not singular, then the 

constraints are second class |34l] . 
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A B 
C D 



In our case, the matrix of Poisson brackets can be written as a supermatrix of the form 
where A is an 8 x 8, -B an 8 x 6, C a 6 x 8, and £) a 6 x 6 matrix. Matrix A has the form 

^ = ^(0)^^(1) ^ 



(19) 



(20) 



where A^^^ contains only bosonic fields and A^^' is proportional to p. The explicit form of the matrices above is given 
in Appendix [XI In order to see whether the matrix is singular, we have to consider the superdeterminant [40l. l4l| 



sdetj^^'^), I = deiA[dei {D-CA-^B)\ 



Here we have 



and 



AetA = - (1 + 2p) 



det {D - CA-^B) = 1 



(21) 



(22) 



(23) 



such that the matrix of constraints is not singular. Hence we have only second class constraints. 

Since {(^(°\ i;^^'')} is not singular, it is possible to obtain the Lagrange multiphers from eq. (fTT]) by considering its 
inverse 



Then, the equation of motion for an observable can be written in terms of the Dirac bracket [Sj 



where 



(24) 



(25) 



(26) 



In order to obtain the inverse of the supermatrix of the constraints we determine first the inverses of A^^-* and £), 
and form the following matrix 



AT 



D 

Multiplying the matrix of the constraints ((T9)) by A/", we have 

y,^''\(f)^^^^N ^ 1-R 

where 



R 



C A 



1(0) 







(27) 



(28) 



(29) 



Then, since R contains Grassmann fields, the inverse of the matrix of constraints is achieved with a finite number of 
powers of R: 



= N(l + R + R^ + R^ + R'^ 



(30) 
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Having the inverse of the matrix of constraints, it can be readily seen that the commutations relations ([4]) are 



reproduced by the Dirac brackets: 



ilx^'^^Xj*! . Hence, the Lagrangian (|12p and the set of 



constraints (|M)) lead to the algebra ([5]). 

The path integral quantization can be performed starting with the Hamiltonian, and integrating over the canonical 
fields and the corresponding momenta. As shown by Senjanovic [13, 113], in the case of second class constraints, the 
path integral is as follows. 



Z = j VXVn\{5[^('^^] sdeti/2{0W,0W| e'/dt[nx-//] 



(31) 



In the present case it is possible to integrate over the momenta, leading to 



Z 



VX 



a(2) 



a(9) 



(l + 2p,)^ e- 



(32) 



where the action in imaginary time is given by 



S 



E 



{l + p,)u,-l 



X7^x; 



\ (xfX^^ + X^'^Xf) +H{X)\ 

i,cr ) 



(33) 



with [3 = l/fc^T the inverse temperature. The term coming from the superdeterminant gives just a shift of the 
chemical potential, and can be ignored. Since X^^ does not enter in the action, we can integrate over it, eliminating 
the constraint 

As a final step, we perform the change of variables introduced in (30l.[33l| that will be useful for further considerations: 



(34) 



X++ = (l-p)(l + f},)/2, X— = (1 - p) (1 - f7,) /2 , 

x+- {1 ~ p) {a, - tcty) /2 , X-+ = {1 - p) {n, + lily) /2 , 

x°+ = VI , = C , 



where we introduced Grassmann field and ip^. Accordingly, we have p — '0+V'+ + V'lV'-- Since we have only 7 
new variables, chosen in such a way that (p^^^ is automatically satisfied, we integrate over X taking care of the 
constraint (f)^-^'^ before performing the change of variables. After the change of variables, the remaining constraints 
look as follows: 



(10) 



For (/)(2) the factor (1 - pf can be absorbed in the chemical potential, such that it reduces to 
After making the change of variables, we finally have 

rl3 



S = 



^ ii-p,)ii-p,){n.-n,-i)-pJ2pA , 

<i,j> i ) 



(35) 

(36) 
(37) 



(38) 



with constraints 4>^^\ (|36p . and (|37p . For the undoped case, the action above reduces to the one corresponding to a 
quantum Heisenberg antifcrromagnet, as obtained e.g. using coherent states [39| . 
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III. STAGGERED CP^ REPRESENTATION 



Here we concentrate us on the limit of low doping, in order to study the consequences of doping on the antifcrro- 
magnetic state present in the undoped case. In such a situation, we can assume a large correlation length for spins, 
such that a long-wavelength expansion is justified. In this limit we can also neglect terms ^ piPj obtaining thus a 
bilinear form in the fermionic degrees of freedom. Then, the action has the following form 

S = Ss + Sf, (39) 

where 



Ss = j '^ dr<j - ^ ^ A[^l,\ ■ drn, + ^ ^ • I , (40) 
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<jj> 



is the action of a pure Heisenberg model. A is the vector potential of a magnetic monopole, (V x yl) • = 1, where 
derivatives are taken in the fi-space. In this case it is given by 



(41) 



(1 + r!,)' (1 + r!,)' 

The fermionic part is given by 

^ \ i.a <i.3> <<i,j>> 

+t" ^-'^.'^ - 1 E ^lA.^^ • + M E 1 ' (42) 

<<<i,j>>> <i,j> i^a ) 

where we allowed hopping to 2nd. nearest neighbors along the diagonals of the square lattice, denoted by << i,j >>, 
with hopping amplitude t' , and to 2nd. nearest neighbors along the principal axis, denoted by <<< i,j >>>, with 
hopping amplitude t" . Such hopping terms have been taken into account, since there is consensus that they are 
present in the real materials [43]. Furthermore, we have to take into account the constraints ([55)) and (P7)) . 

A. Rotating reference frame and staggered CP^ 
representation 

In order to take into account the constraints (|36p and ([37]), and in order to work with smoothly varying fields, we 
consider now the fact that the field Jl is staggered on nearest neighbors and define a local quantization axis for the 
fermions with a rotation that fulfills 

Uln,-aU, = i-iya\ (43) 
where cr", with a = x, y, or z are the Pauli matrices. For i even the condition above is accomplished hy U G SU(2), 

where z z = 1 , and 

n^^za^z. (45) 
For i odd, on the other hand, we take rotation matrices (|44p and 

n'^ = z^al^a^^^ay,z; . (46) 

This ensures that is fulfilled. 
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1. Constraints in the rotating reference frame 



Here we discuss the transformation of the constraints ((36|) and (pT]) on going to the rotating reference frame 
introduced above. For simpUcity of notation, we redefine them as follows. 

^^^^ ^ -^+TTT7f + ^-' 

</5F,2 = -i^+{n^ + my) + i)^{i + n,) , (47) 

Since (pp2 is a fermionic constraint, we can use a corresponding (5- function for Grassmann variables ^ and 

<5(e-n = -(e-?') , (48) 

such that 

S {^*Fa) = {^*F,2) . (49) 

where we defined <fp 2 = (^^ ~ if^j,) + "(/'* {1 + ilz)- Furthermore, going back to the original formulation in terms 
of the X-fields, we can rewrite 0'^^^ in as follows. Insertion of (/)(^^ into cf)^^") leads to 



such that 



and 0*^^^ can be brought to the form 



= X+-X-+ - X++X— , (50) 



X- = ^^£1: , ,51) 



0(9') ^ , (52) 

that after the change of variables we can write as 

<^F,l = (1 - f^.) - V-* (f^x + l^y) ■ (53) 

Then, we can consider the first constraint in (I47p as part of 

ip*p {1-n- cr) = . (54) 

We can introduce now new fermions x — U'^ip and x* — '4'*U, such that 

>f*pU = x*U^ (1 - • <t) [/ = X* (1 T fT") = , (55) 

where the upper (lower) sign is for even (odd) sites. 

In the same way, we have for the second constraint in (|47p 

U^'ipF = U^l-n-(T)Ux^ {lTcrlx = ■ (56) 
However, since the constraints are actually given solely in terms of (f>*p and (pp, we introduce 

e* = ^*pU , e = jjVf , (57) 

such that for the original form of the constraints we have 

^ ■5{eiz; + 9;z,)5{z2e, + zie,) . (58) 



We consider now the action of the constraints on even and on odd sites. 
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i) Even sites. 

ei = d,=0, 9; =2x1, e^ = 2x_. (59) 

Furthermore, we have 

l + n, = 2\z,W (60) 

such that finaUy, 

6 5(^^.2) = s {e*z,) 6 {zie,) = 26 (xl) S (x_) . (61) 

i) Odd sites. 

01 = 2x1, 0,=2x+, 9;=9., = 0. (62) 

In this case we have due to 

l + n, = 2\z2W (63) 

such that finaUy, 

S {ip*p^2) 6[^p^^) = S {9lz;) S {z^9,) = 25 {x*+) S (x+) . (64) 

The constraints above lead to X- = on even sites, whereas x+ = on odd sites. We can therefore work with 
spinless fermions XA on even sites and xb on odd sites, where A and B denote the two sublattices, such that the 
constraints on fermions are exactly taken into account. 

2. Slow CP^ variables in the rotating referenee frame 
After introducing cells j , each one containing one even (A) and one odd (B) site, we can define new fields 

5. = I (-f + 4) ' < = I (-f - 4) ' (65) 

where a is the original lattice constant. Due to the constraints | p=| z^ p= 1, the new fields are subjected to the 
constraint 

zC + Cz = 0. (66) 

This condition can be used to fix the phase of C with respect to that of z, such that they change by the same amount 
under a gauge transformation. This will be discussed in more detail in Sec. IVI Al Equation (|66p also implies that 
both fields are subjected to the constraint 

zS + a^ CC=^- (67) 

We then introduce new fields 



z = z^l-aKC, (68) 

such that the constraint (l67|) is satisfied with zz = 1, and the constraint ([66|) translates into 

zC + Cz = 0. (69) 

The fact that z^ — z^ is of 0{a) can also be seen by going back to the vector representation, where we have 
~ zcTz while Jl"^ + Jl^ ~ C,(tz + ztjQ, i.e. the field Q is directly related to ferromagnetic fluctuations within 
the unit cell. 
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IV. LONG- WAVELENGTH EXPANSION 

Once we have identified smoothly varying fields and their slow modes, we perform an expansion of the action in 
powers of the lattice constant a up to second order, after a transformation to the rotating reference frame. For clarity 
of the presentation we deal first with the action 3$ in eq. (j40p and then with Sp in eq. 



A. Spin action in the staggered CP^ representation 

Here we consider the pure Heisenberg model as given by ([^0]) . We pass to CP^ variables using (1^5]) for even sites 
and (|46)) for odd sites. 

For the Berry phase we have 

- ■ drn, - 2a ^ [zj drQ + drz,] + O (a^) . (70) 

* i 

Terms containing a total time derivative were discarded due to periodic boundary conditions in imaginary time. 

For the interaction term we first discuss our convention in defining new coordinates for the units cells containing 
sites A and B. On passing to the new coordinate system we choose 

^ -^{x + y) , y' = ^{y~x) , (71) 

such that in the new coordinate system the basis vectors for sublattice A and B are, respectively, 

XA - (0, 0) , XB= (a/V2, -a/V2) . (72) 



Introducing the notation 

Gj EE 2iZjaal^Cj(3 , = iZj^alpd^Zjp , (73) 

where fi = x or y, we obtain after a lengthy but straightforward calculation, 

^ . = 8 / da;2 1 2 [G,G* + F.yF;^ + F^^F*^ - {F^^F;^ + F^^F*^ 

+^ [iF,y - F,^) G* + {f;^ - F*:)] | , (74) 

where constants terms were discarded. 

B. Fermionic part in the staggered CP^ representation 

Here we consider the action in the rotating reference frame, i.e. with fermions as defined in Sec. IIII A II 

X = UH, X*=rU. (75) 
After applying the constraints (|36p and (I37p . we can define a new spinor per unit cell 

following the discussion in Sec. IIII A II 
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1. Temporal derivatives 



Defining 



we have 



(77) 
(78) 



where we introduced 



(79) 



At this point wo can check gauge invariance. Since XA and xb have opposite charges, a gauge transformation leads 
to 



XA ^ e+^'^XA , XB ^ e-^^XB 



(80) 



It can be easily seen that the term with the covariant derivative is invariant under the gauge transformations above, 
and that is gauge invariant due to the constraint 



C„ 



-zd^C + Cd^,z + d^cj,{-zQ + Cz) 



(81) 



2. Kinetic term 

For the kinetic term of the fermions we have the contributions proportional to i, t', and t" , that we treat separately 
in the following. We consider first contributions from nearest neighbor hopping, where, after Fourier transformation, 
we have 



(82) 



where 



S(fci, fc2) 



(83) 



with cr± = (cr"^ ± icjy ) /2. The explicit expressions for sW(fci,fc2) are given in Appendix [B) They contain only 
contributions up to ©(a), since, as seen above, terms proportional to t are off-diagonal, and hence, no contribution 
linear in S appears in the final action, as shown in Sec. [V] They enter in a quadratic form after integrating fermions 
out. 

Next we consider contributions from second nearest neighbor hopping along the diagonal. In this case we have 



(84) 



where 



7+ + 



with 



ei (fc) = — 2t' cos {^/2kxaj + cos {\/2kyaj 



(85) 



(86) 
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and 7^ = (1 ± cr^) /2. Again, explicit expressions for ^^^\k, k') are given in Appendix [Bl Contrary to the previous 
case, we consider here terms up to O (a^), since terms proportional to t' connect sites within a sublattice, giving rise 
to diagonal contributions to the self-energy of the fermions, as seen above. 

Finally, we have the contributions from second nearest neighbor hopping along the principal axis: 



where 



with 



-t" ^-^fc- ^ E xl.^k.M) 



Xk2 



<<<i,fc>>> 



ki ,^2 



$(fc,fc') = 62 (k) Sf^ f^> 1 



(87) 



£2 (fc) = — 4t" cos yV2kxaj cos yV2kyaj , 
and k') as given in Appendix [Bl Here again, contributions up to O (a^) have to be taken into account. 



(89) 



3. Spin interaction dressed with fermions 



For such contributions we have 



ki ,fc2 



where 



T(ki,k2) = T+l+T-(7^ 



with T'^ and T displayed in Appendix [B) 



V. INTEGRATION OF FERMIONIC DEGREES OF FREEDOM 



(90) 



(91) 



Once the gradient expansion was performed, we integrate out in a first step the fermionic degrees of freedom, that 
are considered at all wavelengths. In a second step we integrate out the magnetic fast degrees of freedom in order to 
obtain the effective action for the slow magnetic modes. 

Collecting all results obtained in Sec. lIVBl we have the following form for the action 



Xk' 



k,k' 



where k = {ivmk), being Matsubara frequencies for fermions, and we defined 

Gq i(fc, k') = i {lUr, - [e (fc) + ,7 + /i]} 5kk' , 
i.e. the free propagator for fermions in the present theory. The self-energy S is given by 

S(fc,fc') = i(7'-Ar{k,k') - aCr{k,k')+tE{k,k')+t'^{k,k')+t"^{k,k') - jT{k,k') . 
The dispersion relation for the holes is given by the ones in ([55)1 and ([55)1 . 

e (fe) = ei (fc) + 62 (fc) — —2t' cos (^V2kxaj + cos (^V2kyaj — At" cos {^/2kxaj cos a 



(92) 



(93) 



(94) 



(95) 
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This dispersion is obtained for the magnetic Brillouin zone. By rotating to the original Brillouin zone, we have 



and the dispersion is given by 

e{k') 



— ^ i^X ky) , ky — ^ {kx + ky) , 



-At' COS {kxo) cos ikyo) — 2t" [cos {2kxa) + cos (2/cya)] 



(96) 



(97) 



It should be noticed, that the free dispersion of the dopant holes is determined here entirely by t' and t" . This 
is to be contrasted with calculations of the single hole dispersion of the pure t-J model in the self-consistent Born 
approximation [3, |4^, that are in very good agreement with quantum Monte Carlo simulations [l^], where a 
dispersive band is found also in the case t' = t" = 0. Instead, in our case, on passing to the continuum limit, such 
a dispersion can be obtained only by explicitely introducing finite values for those hopping amplitudes. In fact, the 
dispersion of the single hole found previously can be reproduced by ([57)1 for appropriate values of t' and t" {t' = — 0.3t 
and t" = — 0.115i for J = OAt in Ref. [i^l). Hence, on focusing on lower energy scales, operators corresponding to t' 
and t" should be expected to be generated even if such operators are initially missing. 
Arranging terms according to powers of the UV cutoff a, we have 



I](fc,fc') = aI](i)(fc,/fc')+a^S(2)(fc,fc')+C'(a3) , 



(98) 



such that after integrating out fermions and keeping contributions up to O (a^), the fermionic part of the action goes 
over into 



Sf 



-Tr(GoS + iGoSGoS) + 0(a3) 



The terms entering ([M)) can be regrouped as 

with i — 1,2. This leads for the terms linear in E to 



TrGnE 



TrGn 



aE 



(1) 



E« a' 



2v(2) 



a-'E 



E«a- 



and for the terms quadratic in E to 



TrGoEGoE = a^Tr 



GoE^'^GoE^'' + GoE(i)GoE(i) + 2GoEV'^GoEL'^ + O (a^) 



(99) 



(100) 



(101) 



(102) 



Introducing indices t, i, t' , and t" for the temporal, and different hopping processes, respectively, we arrange the 
different contributions in ([M]) as follows: 



^0,1 (^'^') — ^0 



(103) 



while E^-* contains only contributions proportional to t. 

An explicit evaluation using the results in Appendix [B] shows that there are no contributions in first order in a. 
The contributions in O (a^) coming from TrGoE are as follows: 



TrGoS(2'^' 
TrGoE(2'*') 
TrGoS(2'*") 



zdrC + Cdrz) , 

16t"pi J dr J d^x {dxzdxZ + dyzdyz) , 
IQt" p2 J dr J d^x {dxzdxZ + dyzdyz) , 

TrGoE(2.^) = -SJp J dr J d^x i^2GG* + 2FyF* + 2F^F* - {F^F* + F,^F*) 



+V2 



{Fy - F^) G*+G {f; 



F*) 



(104) 
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where p is the density of holes, i.e. 



pi is defined as 



and 



P2 



pi = ^ n (fc) cos (^\/2kxa^ 



— n (k) cos (^V^kxO^ cos a 
k 



In the expressions above, and in the following, n (k) is the Fermi distribution function 

1 



n (k) 



exp{/3[e (fc) -/.]} + ! 



(105) 



(106) 



(107) 



(108) 



We focus now on the contributions quadratic in E, given in eq. (jl02p . The different terms from the expressions 



containing Eq^^ are: 



^0 + ^0 



(109) 



The same expression results for the terms containing E 



(1) 



In the following we consider the different nonvanishing contributions. From the temporal part we have 
TrGoE(i^^)GoE[,i'^)+TrGoE(i'-)GoE(i'-) = ^ dr [ d^x{zdrzf , 



where the electronic compressibility is defined as usually, 

ldp_ 
~ pdn' 

From terms proportional to t we have 

Tr GoE^^^'*^GoEL^'*' = 16fp f At [ d^x [2kiG*G + V2niG* {Fy - F^) + V2ki {F* - F^) G 



(110) 



(111) 



where we defined 



pd^iN 



k 



2 V2. \ JV2\ 
cos I —k^a 1 cos I —kya I , 



and 



K2 



(112) 



(113) 



(114) 



The terms proportional to t' lead to 



TrGoE[,^'*''GoE[,^'*'^ +TrGoE(i'*')GoEii'*') = -32t''^a'^pK3^ j dr \(zd^zf + {zdyz) 



(115) 



where we defined 



A; 



Prom the terms proportional to t" we obtain 

TVGoS(^'*")GoS[,^'*")+TrGoE(i'*")GoSii'*") = -QW^'^-pK^y^ t 

with 



dr 



1 d 1 

K4 = — 



Finally, there is a contribution proportional to t't" of the form 

TrGo4^^* ^Go^^'* ^ +TrGoE(i'*')GoS^i'*") = -Q^t't" a"^ pnsj^ dr (zc^^z) V (^9^^ 

where 

= ^^^^^'^(fe) sin^ ^v^/c^a^ cos ^V^/cytt^ + sin^ ^V^fcyC^ COS ^V^/ca^a^ 
A; 

Collecting all contributions after integrating fermions out, we have 



where 



I 25 K 

dr / d^a; < 16 {t' pi + t" P2) {dxzdxZ + dyzdyz) — (zdrz) 



+16p {J - t'K2) {F*F^ + F;F^) - 8p ( J - 2t2«i) (F^Fy + F^F^) 



contains only contributions with the 2:-field, and 



.(C) 



/3 



dr / d^X J - ^ {zdrC + Cdrz) 



+8p (J - 2fti,) {2G*G + V2 [G* {Fy - F,) + (f; - f;) G] } | , 



contains contributions with ^-fields. 
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(116) 



(117) 



(118) 



(119) 



(120) 



(121) 



(122) 



(123) 



VI. INTEGRATION OF FAST MAGNETIC MODES 

The effective action at this stage contains only the slow modes described by the 2;-fields and the fast modes 
corresponding to the (^-fields, The next step is to integrate out the ^-fields. However, due to the gauge freedom 
introduced by the CP^ fields, a gauge fixing is necessary. 



16 



A. Gauge fixing 



In order to integrate out the C-field we have to discuss the gauge-fixing for the z-fields and its consequences on the 
(■-fields. As already mentioned in Sec. IIII A21 the condition can be used to fix a global phase of the (-held with 
respect to the one of the z-field. In order to see this, we can use the following parametrization of the 2;-field: 



cos (I) exp{-^(f-A)} 

sin(f) exp{-|-z(f +A)} 
where < 9 < For ( we use a similar parametrization 

' pi exp{-i(|-r)} 

c = 

_P2 cxp{+i (i + r)} 

Then, the condition (j66p can be fulfilled with F = A + (2m + l)7r/2, with m integer, that is equivalent to 



(124) 



Ci C2 



Ci I IC2 



I ^1 I \Z2\ 



that can be translated into 



I ^1 II Ci I 

I Z2 II C2 I 



Then, the condition (|66p can be enforced as follows: 



6 (zC + C^) 



^2||C2 



Zl I iQz2 + 4C2) 



Z2 II C2 
Ull 



Cl 



where we solved the constraint (|127p in favor of | (i |. 

Once we enforced the constraint (|66p . we impose the following gauge fixing for the z field 



z-i + z-i 







(125) 



(126) 



(127) 



(128) 



(129) 



as normally done in a CP^ theory. The implementation of the constraint (|128p in the part of the action containing 
(-fields (|123p is given in Appendix O 



B. Ciiange of measure 



We discuss here the change of variables from the real vectors fl^'^ to the complex fields z and (, taking into 
account the Jacobian of the transformation, together with the constraints for the z- and (-fields already discussed in 
Sec. IVI Al For this purpose we first consider and 



Then, we have 



ZjJl- a^CjCj - aQ ](t[zjJ1- a^QQ - a(. 



3 ' 



z-j y 1 — ci?CjCj + a(j (T^ (T z* W 1 — a'^QQ + aC 



(130) 
(131) 



Therefore, 



A |2 



B |2 



ZjZj (1 - a^CjCj) + a^QQ - a {z^Q + Qzj) J 1 - a^QC 



ZjZj (1 - a^CjCj) + a^QCj + a {zjQ + CjZ]) y/ 1 - a^CjCj 



6{\n^\ -1) 5 (1 I -1) ^ (^1 + ^a^CC^ S {zz - 1) S {Cz + -zC) 



(132) 



(133) 
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The transformation of the last constraint due to gauge fixing was already discussed in (I128p . On the other hand, we 
have due to p?7)) 



cc 



1^1 1= 



(134) 



where we used zz = \. 

Due to gauge fixing, and without imposing the constraints on the moduli of the vector fields $7"^ and Jl^, since 
we take them into account with (|133p . we have independent variables Imzi, Rez2, IniZ2, | Ci li R'6C2j and ImC2. The 
Jacobian going from vector fields to CP^ variables is in an expansion in powers of a, 



where 



and 



Jo = 32 



J^a" (Jo + a' J2) + O {a^) , 
zlmzi(2|zi |2-1) ((2*^2 + ^(2) 



zi 



J2 = -48 



IIC2 I ' 

zlm^i (2 I zi |2 -1) (C*Z2 + ZIC2) I C2 



Zl 



Z2 



such that we finally have 



^ _ 3oo (2 I ^1 P (C2*^2 + ^2X2) A 32IC2I 

J ~ a 6Z I I I 11,1 I -L I r, |2 



Z2 II C2 



Zl 



(135) 



(136) 



(137) 



(138) 



Taking into account the Jacobian above, the constraints in Sec. IVI A[ the transformation of the constraints (|133p 
together with (|134p . the measure changes as follows 

j dn^ 5{\n^\ -1) <5 (1 I -1) ^ j dimzi dRez2 cllmz2 5 {zz - 1) '^""^^ ^^^J 

X / dReC2 dImC2 ■ (139) 



C. Effective field theory for the magnetic properties of the t-J model 

The final step is to integrate out the C-field. Therefore, we concentrate on the corresponding integrals. 



I?ReC2 I?ImC2 e" 



?(C) 



j PReC2 CTmC2 ^xp j - dr j d^x + A*C| + AC2*' + (S - T) C2* - (S* + T*) C2] | , 



where we introduced the following notation: 



Zl 



A = 8Jz. 



2 J 

2 (1 - 2~p) 



Z2 



Zl 



Zl drZi 



Z2 



2V2J 



Z2 



^1^2 
Z2 P 



{Fy - F.,) + {f; - f; 



T2 ^2 dTZ2 



Zl P 



(140) 



(141) 
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We go now over to real and imaginary parts (2 = V + such that for eq. (|140p we have 







j Vrj cxp 









(142) 



with = A + A* + A, 6 = A - (A* + A), C = A* - A, X> = S - S* - (r* + r), and £ = - (S + S*) + (r - T* 
Performing the integrals over the C-fields leads to the effective action 



(143) 



Next we proceed to evaluate the different contributions to the action. We consider first the term that affects the 
measure 



AB + C^ = [j + 8p(j-2tV)] (1-2 I zi 



Then, 



exp 



-^\n{AB + C^ 



zi 



[j + 8p{j-2t\i)] (1-2 I zi p)| 



such that the measure and constraints of the path-integral are now 

Z I VzVz5{zz-l) ^(zi + z*) (z^-z*) 



that are the ones usually appearing in the CP^ model. 
After a lengthy but straightforward calculation we have 

[AE"^ - 2CV£ - BV^) 



(144) 



(145) 



(146) 



4(y^e + C2) 



{2J 


[(1 


-2/5)1 


2 






a 




1 


[(1-2/5)] 




a 




-J 


dxzdxZ 4 


dy 



drZ drZ + {zdrZ) 

[drZ {dyz - dxz) - drZ (dyZ - dxz)] 



-dxzdyZ — dyzdxZ — 2 {zdxz) {zdyz) 



(147) 



Taking into account the contribution to the action in (|122p . that remains unaffected by the integration over the 
C-fields, we arrive at the effective action for the z-fields 



with fi = T, X, y, and 



S = /drd^x^— df_,zdf,z + jf^{zd^zy 



2Jo2 



(148) 



(l-2p)' ' 



gx=9y = [j(l + 8/5) + 16/5i2(Ki-K2) + 16(t'pi+t"p2)] ' 
4pKj 



Ir = 1 



(1-2/5)^ 



J(l 



16p 



(k1 - K2) + t K3 + 2t" Ki + At't"K5 



[J (1 + 8p) + 16pt2 (ki - K2) + 16 (t'pi + t"p2)] 
where J was defined in (|14ip . Equation p48p together with p49p are the main result of this work. 



(149) 
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VII. DISCUSSION OF THE RESULTS AND CONCLUSIONS 



The action (I148 [) at which we arrived, is of the form previously analyzed in the context of frustrated quantum 
antiferromagnets [l2l [l^ . In the absence of doping, i.e. setting p — pi — p2 — have — Ix — ly = such 

that the model reduces to the CP^ model appropriate for a Heisenberg antiferromagnet, making thus explicitely the 
connection with the undoped case. As is well known, the excitations in that case correspond to bosons describing the 
transverse fluctuations of a vector field pertaining to the 0(3) non-linear cr-model. From the point of view of the CP^ 
model this corresponds to a phase where the z-fields are confined. 

However, in the presence of doping, when the couplings 7^ 7^ 1, the model does not correspond any more to a 
collinear antiferromagnet but describes in general coplanar incommensurate quantum antiferromagnets (l^ . [l3| . This 
can be seen by constructing R e SO (3) out of matrices g e SU(2), as follows: 

Rab[g) = iTraV's^ , (150) 

where g is given in terms of the z-fields like in (|44p . Then, it can be readily shown that the action (|148p can be 
expressed in terms of SO (3) fields: 

S ^ ( AtAx'^Y -^Tr {d^Rd^R-^ - ^^d^RQD^R-^) , (151) 
J ^ 8g^ 

where Q — diag(l, 1, — 1). Expressing the matrix R as R — (rii, n2, ''t-s), where the vectors rii fulfill nf = 1 for 
i = 1,2,3 and rii ■ rij = for i j, the action psip takes the form 



S = I drdx^Y — 



2 



Y ^fj^ni ■ df^rii - (1 + 7p) (ni • 5pn2)^ 



(152) 



showing that in general a coplanar configuration is favored. In the case 7^ = 1, the action above reduces as expected 
to the 0(3) non- linear cr-model by virtue of the relation 

a^na • a^ng = a^rxi • S^ni + a^ria • ^^"2 - 2 (ni • 9^n2)^ • (153) 

On the other hand, for the general case 7^ 7^ 1, assuming for simplicity that 7^ = 7 and g^ — g for fi — t,x, y, it 
can be seen that, following Ref. [l3|, the gauge fields responsible for confinement in the case 7 = 1, acquire a mass 

M^ = -^(l-l). (154) 

Moreover, it is expected that in the infrared limit, the behavior is dominated by an 0(4) fixed point [l^, [l^li 
it is expected that the coupling 7 scales to zero, raising the mass of the gauge fields. Further arguments in Cavour 
of deconfinement were advanced by showing that for 7 7^ 1, the mass of the gauge field is inversely proportional 
to the square root of the temperature jl4| . and therefore, a confinement-deconfinement transition as a function of 
temperature can be expected when 7 departs from 1. Hence, we see that the possibility of deconfined (bosonic) 
spinous is opened due to doping of a collinear antiferromagnet. 

Finally, we would like to remark, that the present results show some differences from those obtained previously 
by one of the authors [48] and collaborators, dealing with a doped antiferromagnet modeled by the so-called spin- 
fermion model. This model is the starting point that leads to the t-J model in the limit where the exchange coupling 
between the dopant hole, that mainly resides on the oxygen orbitals, and the copper hole form the Zhang-Rice singlet 
[27! . The gradient expansion for the spin-fermion model based on the assumption of a short-range antiferromagnetic 
order led to an 0(3) non- linear cr-model, that as a function of doping had a transition to the corresponding quantum 
disordered phase [49]. Hence, under the same assumption as in the present work, no hint to deconfinement of spinous 
was obtained. An SO (3) non- linear cr-model can be obtained from the spin-fermion model only under the assumption 
of an incommensurate coplanar short-range order in the microscopic model j50| . Therefore, although considering 
the symmetries present in a doped antiferromagnet, both models appear equivalent, we conclude on the basis of the 
present results that imposing the constraint against double occupancy leads to a richer picture of the possible phases 
of a doped antiferromagnet, with a minimal number of assumptions. 
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APPENDIX A: SET OF CONSTRAINTS FOR THE t-J MODEL 



We give here the expHcit expressions of the constraints arising from considering the canonical momenta of the 
system, and the matrix of constraints. First we recall the definitions of right and left derivatives in order to reach a 
self-contained presentation. 

Given generators of a Berezin algebra with a = 1, . . . , fc, the right derivative is defined as 



dr 



i=l 



(Al) 



where P^a) =0,1 depending on the parity of the generator (even or odd, respectively). For later use, we introduce 
also left derivatives defined as 



A. 



k 

E 

i=l 



(A2) 



In both cases it is understood that when the index j runs beyond the interval [1, fc], then P(aj) = 0. Both derivatives 
are equivalent to ordinary derivatives when dealing with even generators z". 
With the rules above we obtain the following constraints from ([T5)) . 



0(3) = n++ , 
= n— , 



(b^^^ = n+- + 
= noo , 

(13) ^ jjO+ _ 



(1+p)m-1 



2X 



-+ 



^(11) 
^(14) 



n 

^^ 



(1+p)m-1 
^— 4p— 



(A3) 



The matrix of constraints (I19|) has the following components 





/ 





1 



-1 X— 
-X- 



-1 X++ 



1 











V 






































^{2+X + +) 

X + + X + - 








j(2+X + + )p 


X + + X- + 




















Vo 









1 

-X- 



i 

~ 4X + - 
i 







X-+ 

i 

4X+- 





t{i+x-- 

4X++X + ' 








x+- 

i 

4X- + 



~XTT 



i{i+x~-) 

' 4X + + X- + 









i{2+X + + ) 

X + + X + - 



i{2+X--) 
4X+ + X + - 






i(2+X + + )p 

X+ + X- + 



~XTT 



-t(2+X--)p 
4X+ + X~ + 





1 







1-X++ 

















t{i+x--) 

4X+ + X + - 








i[i+x--) 

























1 





-1 


0^ 








i(2+X~~)p 

~ 4X+ + X+- 
i[2+X )p 

4X + + X- + 







0\ 











0/ 



(A4) 



(A5) 



The matrix B contains only single Grassmann generators. 



B 



/ 








x°+ 


x^- 




-X-" 


\ 










x^- 


_X+o 


-X-" 
























B41 





B43 


B44 


B45 


B46 









_X+o 


B53 


B54 


B55 


B56 


































-x°+ 






x-° 




V 




















/ 



(A6) 
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where 



-B31 — 

-B43 — 

-B45 — 

B53 — 

B55 — 



x"+x+- 

(X++)^ ' 
i(l+X--+2X+ + X--)X°+ 

4X++X+- 
i[l+X--+2X+ + X--)X+° 

~ 4X+ + X + - ^ ' 

i[l+X--+2X+ + X--)X°+ 

4X+ + X-+ ' 
i[l+X~'+2X+ + X^^)X+° 

4X++X- + 



B41 

-B44 
-B46 

-B56 



X" 



X++ ' 

i(l+X--+2X+ + X~-)X° 

4X + + X + - 
i[l+X--+2X++X--)X~° 



4X+ + X--X+- 
i(l+X--+2Jf++X-- 



)X° 



4X++X- + 
i(l+X — +2X+ + X" 



4X + + X- + 



-)x-° 



(A7) 



Furthermore, C = B , and 



/ 














x+- 

X + + 


-l\ 












-X++ 





















— i 










-X++ 











—i 




x+- 

X + + 





—i 











V 


-1 








— i 





J 



(A8) 



APPENDIX B: FERMIONIC TERMS IN THE STAGGERED CP^ REPRESENTATION 



We give here the expHcit expressions for the fermionic contributions up to O (a^)- Defining vectors 

= (0,0), = (0,V2a), = V2a(-l,l), a;^^) = V2a(-1,0), (Bl) 

we have for the fermionic terms si*'' (fci, fc2) with i = 1, . . . , 4 

S^^(fci,fc2) = -^exp|i ki ■ xa- k2 ■ (xb + x^'^^^^^exp[i{ki - k2) ■ Xj]Ijlj'\g , (B2) 

j 

where Uj^\g denotes the matrix element of the product of the SU(2) matrices in nearest neighbor sites, connecting 
sublattices A and B within the unit ccU j. The index i denotes the four nearest neighbors starting from sublattice 
A. Using the notation introduced in ([75]) . the gradient expansion of the products of the SU(2) matrices up to 0{a) 
have the following form 



(B3) 



For the fermionic terms ^^^(fci, k2) we have 



S^yfci,fc2) = exp{i [fci • a;s - fc2 • - a;0]}X!*^^P ^ ^2) • a^j] ' (^4) 



with U% = U%. 



Next we list the terms originating from contributions proportional to t' . For ^^^Cs^i, ^2), ^ = 1, . . . , 4, we have 

'^^AA iki,k2) = — exp [±ik2 ■ X4] exp [i (fci - fea) • a;^] J/j^^/j , (B5) 



^aA"' (fci : ^2) = exp [Tife2 • 3^2] (^1 ^ ^^2) • Xj] U^yi , (B6) 

j 

where Uj^\j^ denotes the product of SU(2) matrices on second nearest neighbor sites on sublattice A along the 
diagonals, starting from the unit cell j, and i numbers the four possibilities. They have the following expansion in 
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powers of a: 



Zjdlzj - \/2 {zjdyCj + CjdyZj) 



U^^AA — V2azjdyZj + 
Uj^AA = -V2aZjdxZj+i 



■jdlzj + V2 {zjdxCj + CjdxZj) 



Uj%A = -"^azjdyZj + [zjd^Zj + V2 {zjdyCj + CjdyZj)^ . (B7) 

For ^'3^3, we have the same phase factors as for ^f^^ but we have to replace Uj^aa ^j%B ^^^^ foUows: 

^i,SB = ^/^azjdxZj + zjO^Zj + V2 {zjdxCj + QdxZj) 
^j%B = V^azjdyZj + a'^ zjdpj + V2 {zjdyCj + QdyZj) 



(3) 



j,BB 



-V2azjdc 



Zj + a 



Zjdlzj - y/2 {zjdxQ + QdxZj) 
Zjdpj - V2 {zjdyCj + CjdyZj) . 

For the contributions proportional to t" we have for processes involving the sublattice A, 



^j%B = -V2aZjdyZj 



^AA {^1,^2) = —exp[±ik2-X3]^exp[i{ki-k2)-Xj]Uy 



(1,3) 
AA ' 



^A'A{ki,k2) = —exp[Tik2- {X2- Xi)]^exp[i{ki- k2) ■ Xj]U>'^J^l , 



(B8) 

(B9) 
(BIO) 



where C^j*^^ denotes the product of SU(2) matrices on second nearest neighbor sites on sublattice A along the principal 
axes, starting from the unit cell j, and i numbers the four possibilities. They have the following expansion in powers 
of a: 

Uj^AA = V2a{zjdxZj - ZjdyZj) 

W [zjdlz-j - 2zjdxyZj + Zjd^Zj - y/2 [zj {dxQ - dyQ) + Q {dxZj - dyZj)] } , 



U. 



(2) 



j,AA ^ y/2a{zjdxZj + ZjdyZj) 

W {zjdlzj + 2zjdxyZj + Zjd^Zj - V2 [zj {dxQ + dyQ) + Cj {dxZj + dyZj)] } 
^■%A = -V2a{zjdxZj-ZjdyZj) 



U 



(4) 



W {zjdlzj - 2zjdxyZj + Zjd^Zj + \/2 [zj {dxQ - dyQ) + Cj {dxZj - dyZj)] } , 

— V^a ( Zj dxZj + Zj dy Zj ) 

W [zjdlzj + 2zjdxyZj + Zjd^Zj + V2 [zj {dxCj + dyCj) + Cj [dxZj + dyZj)] | . 



(Bll) 
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For $^^5, we have the same phase factors as for 3"^^ but we have to replace f^j*!^ by J/j ^^^^ follows: 

Uj^BB V2a{zjdxZj - ZjdyZj) 

W [^zjdlzj - 2zjd^yZj + zjdpj + V2 [zj {d^Q - dyQ + Cj [d^Zj - dyZj)] | , 
f^j.BS = V2a (zj-a^j^ + 

^Zjdlzj + 2zjdxyZj + ZjdyZj + V2 [zj [d^Cj + dyCj) + Cj (5;^% + dyZj)] | , 

s = ^ (^i % ~ ^t/ % ) 

W [^zjdlzj - 2zjd^yZj + Zjdpj ~ V2 [zj {d^Cj - dyC]) + Cj {dxZj - dyZj)] | , 
^j%B = -V2a{zjdxZj + ZjdyZj) 

W [^zjdlzj + 2zjd„jZj + zjdpj - \/2 [z^ [d-Xj + dyQ + Cj (d^Zj + dyZj)] | . (B12) 

Finally, we list below the contributions from spin interactions dressed with fermions. We defined 

^ ^^('='(fcl,fe2)±^('''(fcl,fc2)l , (B13) 



with 



N 

where 



4„2 

^(^■°)(fci,fc2) = —e^p[t{ki-k2)-XA,B]J2exp[i{k,^k2)-x,]T,, (B14) 



T, = 4 {2G,G* + 2^;,^- + 2F^^F;^ - {F^^F*^ + F^^F*^ + V2 [{F^^ - F^^) G* + G^{f;^- F*^)] } , (B15) 
contains the contributions from $1^ • rij. 

APPENDIX C: EFFECTIVE ACTION 5^^' 

We display here, in order to facilitate a reproduction of our results, the different contributions to the part of the 
action containing the ^-fields (|123p after imposing the constraint (|128p . We consider first the part containing temporal 
derivatives in (|123p 

zdrC + C drZ = Z, drZl - zl d^Z^ + Z^ OrC^ - ^2 drQ ■ (CI) 

I Zl |2 I Zl |2 

We consider next the interaction part in (I123p . Here we have 

G = -2zi ( + C2) . (C2) 



Zl P 



Then, for the different terms entering (|123p we have, 



G*G = 4[Lll^^±^\C,\'+zfe,+zlcf 
G*F^ = ^2zl [ ^ + q) F^ , 



Zl 12 



2 /■* 



F^G - -2F:zi [^+(2 ] . (C3) 
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